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Abstract 

We present a detailed, quantitative study of the competition between 
interaction- and disorder-induced effects in electronic systems. For this the 
Anderson-Hubbard model with diagonal disorder is investigated analytically 
and numerically in the limit of infinite spatial dimensions, i.e. within a dy- 
namical mean-field theory, at half filling. Numerical results are obtained for 
three different disorder distributions by employing Quantum Monte Carlo 
techniques which provide an explicit finite-temperature solution of the model 
in this limit. The magnetic phase diagram is constructed from the zeros of the 
inverse averaged staggered susceptibility. We find that at low enough tem- 
peratures and sufficiently strong interaction there always exists a phase with 
antiferromagnetic long-range order. A novel strong coupling anomaly, i.e. an 
increase of the Neel-temperature for increasing disorder, is discovered. An 
explicit explanation is given which shows that in the case of diagonal disorder 
this is a generic effect. - The existence of metal-insulator transitions is stud- 
ied by evaluating the averaged compressibility both in the paramagnetic and 
antiferromagnetic phase. A rich transition scenario, involving metal-insulator 
and magnetic transitions, is found and its dependence on the choice of the 
disorder distribution is discussed. 



I. INTRODUCTION 

The investigation of interacting electronic systems is one of the most intriguing, albeit 
difficult, subjects in condensed matter physics. The same is true for the study of disordered, 
non-interacting electrons. In view of the theoretical complexity of the two problems taken 
separately it is understandable that their combination, i.e. the simultaneous presence of 
randomness and interactions, as found in many real systems (e.g. doped semiconductors 
near the metal-insulator transition, high-T c superconducting materials close to T c , etc.), 
leads to new, fundamental questions to which only few secured answers are known. This is 
all the more true when the interaction and/or the disorder is strong, since there exist hardly 
any tractable, and at the same time controlled, theoretical method of investigation in this 
limit. 

An important starting point for the investigation of interacting, disordered systems was 
the field-theoretic approach developed for the treatment of non-interacting electrons, i.e. the 
scaling theory of Anderson localization (l|-§|. A generalization of this theory to finite in- 
teractions by Finkelshtein ||] provided essential new insight |J. However, the appearance 
of local magnetic moments in the renormalization group treatment discovered by him || 
and Castellani et al. |7[] turned out to be a fundamental obstacle for the study of the metal- 
insulator transition (MIT) itself [|5]||. A microscopic origin of this instability towards the 
formation of localized moments can already be traced within Hartree Fock theory for the 
disordered Hubbard model ( "Anderson-Hubbard model" ) with off-diagonal disorder |To| , pT 



The results indicate that the above renormalization-group approach, as well as is starting 
point, i. e. Anderson-localization, are not appropriate when it comes to the investigation 
of three-dimensional disordered electrons with intermediate or strong interactions close to 
half-filling (n ~ 1). After all, a MIT occurs even without disorder in this case, and the 
lattice periodicity becomes essential. A non-perturbative and qualitatively quite different, 
but essentially uncontrolled, renormalization technique by which disorder and interactions 
can be treated on equal footing is the real space renormalization group approach. It was 



first applied by Ma Jl2| and recently used by Yi et al. [TI| to investigate the Anderson 



Hubbard model with diagonal disorder, i.e. random on-site energies, in dimensions d — 1,3 
at n — 1 and in d = 2 for n ~ 1, respectively. However, in these investigations the formation 
of antiferromagnetic long range order (AFLRO) which will always set in at large enough 
repulsion (at least at n = 1) was not considered at all. The strong coupling limit of this 
model was studied by Zimanyi and Abrahams [jl4| using a slave-boson formulation of the 
corresponding t — J model. The considerably simpler, but still highly non-trivial, case of 
disordered, spinless fermions in d = 1 was also addressed recently [Jl5| , |l6| . 

To obtain a global picture of the properties of interacting, disordered systems it is de- 
sirable to know the solution of a simple, microscopic model which is valid for all input 
parameters (interaction, disorder, temperature, band filling). Since exact solutions are not 
available in d = 2, 3 one would like to construct, at least, a thermodynamically consistent 
mean-field theory that is valid also at strong coupling. Such a (non-perturbative) approxima- 
tion is provided by the exact solution of a model in d — oo. It is now known that even in the 
limit d — > oo [|H],|n| the Hubbard interaction remains dynamical |IU] and leads to a highly 



non-trivial single-site problem P0| , |2"T| | with infinitely many coupled quantum degrees of free- 
dom. This problem is, in fact, equivalent with an Anderson impurity model complemented 
by a self-consistency condition |H],[24],[27]] and is thus amenable to numerical investigations 
within a finite-temperature quantum Monte-Carlo approach [28| . In the absence of disorder 



this technique was already used by several groups to investigate the magnetic phase diagram 
25|j29|| , the Mott-Hubbard transition |30|j32]|29[| , transport properties |H|, and lately also 



superconductivity in a two-band version |34| of the Hubbard model in d = oo. These inves- 



tigations were also extended to the periodic Anderson model [35] and the Holstein model 



The dynamical mean-field theory of interacting electrons obtained in the limit d — oo 
has many appealing features; however, it also has its limitations which are due to the purely 
local nature of the theory. One effect which cannot be described within the d = oo mean- 
field theory of disordered electrons is the phenomenon of Anderson localization. That is, 



the electrons in infinite dimensions are delocalized for arbitrary strength of the disorder. 
The absence of localized states is caused by the above mentioned single-site dependence, 
i.e. by the fact that momentum dependence enters the physical quantities only through the 
dispersion relation whereby vertex corrections to the conductivity vanish identically ||37|| . In 



fact, it was shown P8| , |5§| that for vanishing interaction, when the Anderson-Hubbard model 
reduces to the Anderson disorder model, the well-known and much-used "Coherent Potential 



Approximation (CPA)" [ 40|j41 |, which does not describe Anderson localization, yields the 
exact solution in d — oo. Nevertheless it should be stressed that even in the absence of 
Anderson localization the presence of disorder can have a significant, non-trivial effect on 
the properties of an interacting system. First results of a quantum Monte-Carlo study 
of the competition between disorder and interaction in the Anderson-Hubbard in d = oo 



with diagonal "binary alloy" disorder were recently reported by us |[42| . The importance of 
separable, off-diagonal disorder in the formation of local moments and the Mott transition 
was investigated by Dobrosavljevic and Kotliar |53[ . 



A good impression of how complicated the phase diagram and the dynamical properties 
of interacting electrons with and without disorder are may be inferred from the results for 
spinless electron in d = oo [H1,HH]. The fact that the interactions reduce to their Hartree- 



contribution in this case |19[, i.e. are no longer dynamical, does not at all imply that the 
properties of the model become trivial - quite the contrary. Many details of this model can 
be obtained analytically in d = oo, thereby providing valuable insight into the effects of 



interactions f45| , the interplay between interaction and disorder |43], and the d = oo limit 



itself [46 



The correlations caused by interactions between electrons and by the scattering of elec- 
trons in a disordered medium, respectively, can be very different; indeed, they may lead to 
opposite effects. For example, even an arbitrarily weak repulsive Hubbard interaction be- 
tween electrons on a bipartite lattice will induce AFLRO when the band filling is sufficiently 
close to n = 1 (at least in dimensions d > 2). By contrast, it is the very nature of a random 
potential to oppose order. Hence interactions and disorder are expected to compete with 



each other. This conclusion may appear self-evident and almost trivial but is really only 
based on a simple-minded superposition of the two individual physical effects (which may, 
in fact, be quite correct for weak interactions and disorder). When the interaction and/or 
disorder is strong, however, this picture looses its basis. In general the mutual interplay of 
those two effects can be expected to lead to novel, non-perturbative, quantum-mechanical 
many-body phenomena, e.g. new phases, which have no analog in non-random, interacting 
or disordered, non-interacting systems. 

In view of the interplay between disorder and interaction there are a number of simple 
questions whose answers will be non-trivial, for example 

1. Is there a finite disorder strength above which AFLRO ceases to exist? (After all, the 
fluctuations of the ensembles around the average value n — 1 induced by the disorder 
might become so large that AFLRO is made impossible). 

2. Does disorder always reduce the critical temperature for AFLRO, as one might expect 
intuitively? 

3. Do interactions always drive a system away from metallic behavior into an insulating 

state? 

4. How important is the choice of disorder distribution on the results? 

These are questions which need to be answered by means of controlled investigations of 
a well-defined microscopic model. It is our intention in this paper to provide explicit, 
numerically exact answers to the above questions within the dynamical mean-field theory 
obtained in the limit of high spatial dimensions. A preliminary account was presented in Ref. 
4"2| . The paper is structured as follows. In Sect. 2 we introduce the model and discuss the 



form of the corresponding averaged free energy in d = oo. In Sect. 3 we calculate correlation 
functions, i.e. the averaged compressibility and staggered susceptibility, to determine the 
thermodynamic stability of various phases. The details of the model (input parameters, 
disorder distribution etc.) are specified and the numerical procedure is discussed in Sect. 4. 
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This leads to the construction of the magnetic phase diagram in Sect. 5. In Sect. 6 the 
results on magnetic phase transitions are then complemented by those on metal-insulator 
transitions. A discussion of the combined phase transition scenario closes the presentation. 

II. MODEL AND AVERAGED FREE ENERGY IN D = oo 

The simplest microscopic Hamiltonian for conduction electrons interacting via the local 
Hubbard interaction in a disordered system is the so-called Anderson-Hubbard model, which 
can be written as 

Here cf a (c ia ) create (annihilate) a a-electron on site i, with h ia = cf a c ia and /i CT is the 
chemical potential of the cr-electrons. In general both the hopping amplitude tij and the 
atomic energies ej in (P depend on the configuration of the atoms distributed on the lattice. 
In this paper we restrict ourselves to the case of uncorrelated, diagonal disorder. Hence we 
assume the hopping elements to be independent of the randomness (with tij = —t for nearest 
neighbor sites and t^ = otherwise) while the atomic potentials e, are chosen as random, 
their distribution being purely local; thus short-range order is neglected. The Hamiltonian 
under investigation is then given by 

H = H Hub + J2(ei- fx a )n ia (2) 

where H Hub is the usual Hubbard model 

t 



H 



Hub = — t= Y, d t*Cj* + Uj^nifhil- ( 3 ) 



Here the hopping matrix element was written as [17 



t = t*/VZ , t* = const, (4) 

with Z as the number of nearest-neighbors on the lattice (e.g. Z = 2d on hypercubic lattices 
in d dimensions). This scaling ensures that (^) remains non-trivial even in the limit Z — > oo. 
We can then write the averaged full energy (grand potential) as 



n«, = -/9 _1 <(1b Tr exp(-/?F))> a , (5a) 

where /3 = 1/ksT. Here the average over the disorder is defined as usual as 

(X) av = n / de i P(e i )]X(e 1 , . . . , e L ) (5b) 

where L is the number of lattice sites, and -P(ej) is the distribution of the random potentials 
tj. The explicit form of the function P(ej) will be specified later. 

In realistic dimensions (d = 2, 3) it is impossible to perform the trace over all fermionic 
states in (|2|) exactly. However, in the limit of infinite coordination number or spatial di- 
mensions significant simplifications occur in the evaluation of thermodynamic properties. 
Two of us recently showed |39| that the averaged grand potential of the disordered Hubbard 



model in d = oo can be reduced to an expression where the averaging is performed only on 
a single site. The explicit expression for the grand potential in the paramagnetic phase is 
given by 



" Lav 



/oo 
dEN(E) hi[iu n + \i a - S CTn - E] 

+ /r 1 £inG; n 1 -/r 1 <inZ{s-\ ei }> a „ , (6) 

an 

where N(E) is the density of states of non-interacting electrons, u n = (2n + l)7r/3 _1 are 
Matsubara frequencies and the generalized atomic partition function Z is represented by a 
Grassmann functional integral 

Z{g-\e i } = Jv^V^e A (7a) 

with the local action 



A 



EC„(^ - £i)Vw - U / dr^(r)MrWi(r)Mr) 

a,n Q 



(7b) 



Here ^(t) are Grassmann variables that depend on the imaginary time r, and ^> an are 
their Fourier transforms into tu n -space. Only the random potentials e« are configuration 
dependent. Furthermore, 
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G<rn — G an + Scrn (8) 

is an effective local propagator of the electrons. The grand potential (|6]) is a functional of the 
complex quantities G an and S crn which at this point enter as variational parameters. That 
is, the physical values of G an and £o- n (namely, the local Green function, G an = G« )fm , and 
the self-energy, E CTn = £ji ifm , of the electrons, respectively) correspond to those for which 
(H) is stationary [0,0 . 



Each term of Q av in (|6|) has a clear mean-field interpretation, which allows us to outline 



a simple construction of f2 a „ on a mean- field level p9|JT8[] . We first introduce a homogeneous 



[48|1 , effective, energy-dependent potential S CTn which is defined in such a way that in the 



thermodynamic limit the system of disordered, interacting electrons is equivalently described 
by a system of non-interacting electrons in the potential S CTn . The free energy fi ra {^™} is 
now constructed as follows f|9| : (i) we start with the free energy of non- interacting electrons 
in the potential S CTn (first term in (§)), (ii) then we remove the potential E CTn from site i, i. e. 
subtract its energy contribution (second term in (|6|)), and (iii) replace it by the potential 



Via = \Ufii_ a + 6j — fi a , and finally average over e^ (last term in (|6]) ) [5(J. To determine the 
potential S CTn we demand that the free energy obtained by this construction be stationary 
under variation of Ti an , i. e. dVL av / dY, un = 0. This is a self-consistent equation for E CTn . The 
major advantage of the mean-field grand potential (^ is the property that the variational 
parameters are not explicitly configuration dependent p9 |, i.e. they are determined only 



by averaged quantities which result from the averaging in (||). Since the averaging is local, 
the disorder only leads to local correlations on the same site. Different lattice sites are 
effectively decoupled and all the information about the surrounding sites is contained in the 
effective potential S CTn . Indeed, S CTn is a generalization of the "coherent potential" known 
from the theory of random alloys [[mUD • This may be seen as follows |18| , [51|| : if we neglect 



the Hubbard interaction U, the Grassmann functional integral ( ffal) becomes Gaussian and 
can be performed exactly, leading to 

/oo 
dEN(E)\n[iiu n + fi a -i: (7n -E] 



+ (la{l + G ffn {E m -e i ))) av } (9) 

This is precisely the grand potential of a random alloy with non-interacting electrons within 
CPA, with S CTn as the coherent potential expressed in Matsubara frequencies pl| , |5T| . Hence 
(^) can be understood as a field-theoretic generalization of the CPA to random alloys with 
interacting electrons which becomes exact for Z — > oo. 

It should be noted that the averaged grand potential (§]) does not contain the complete 
information about the equilibrium physics of random systems [[39" |. For example, averaged 



products of Green functions contain new, non-trivial correlations and cannot be derived 
from @ . The averaged grand potential does, however, carry the full information about the 
equilibrium thermodynamics of the model. We concentrate in our paper exclusively on the 
thermodynamics of the Anderson-Hubbard model. 

It is well-known that, for large enough U, the Hubbard model itself leads to AFLRO at 
(or even close to) half-filling. This type of long range order must therefore be considered 
even in the presence of disorder. For the sake of simplicity we investigate only AFLRO on 
bipartite lattices. In this case the breaking of symmetry is caused by a staggered magnetic 
field h s with (+h s ) on A- and (—h s ) on £>-sites. To be able to study this type of order 
quantitatively we extend the averaged grand potential (f|) to the symmetry-broken phase. 
The averaged grand potential of a general antiferromagnetic solution has the form 

R— 1 , poo ,. , 

Q av /L = — i^/ dEN(E) In (iu n + \i Aa - £ A(m )(iuj n + \i Ba - S B(Tn ) - E 2 

-EE ln ^ n + £ (lnZ,) av } (10) 

a,n j=A,B 1=A,B 

where the index 7 = (A, B) = (+1,-1) corresponds to sublattices A, B respectively, and 
/i 7cr = \x + cri^yhs + h), with h as a magnetic field. The local partition function Z^ in (|10"D 
is given by (0) with the replacements ty an — > ^V") ^(t) — > Vv( r )> ^an ~^ G~l n etc., and 
integration over all Grassmann fields, i.e. Z^ — J VipVip* e A ~< with 

A = Y,^*vn(G^ n + S 7CTn - e;)?/vn -U drip^ (r)^ 7T (r)V^ (r)^ 7i (r) . (11) 

" 



rrn 
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The full partition function is given by Z = ZaZb- As described above the physical values 
of G 7(T „ and £ 7<Tn in ( |I0D are then found from the stationarity conditions 

5Q av /5G lan = , 5Q av J '<5£ 7(jn = 0. (12) 

In the following we will set h = in flTU| ) and exclude the existence of /ernmagnetism. In this 
case the average number of electrons on sublattices A and B are equal. Then there appears 
a new symmetry: quantities with indices 7, a and —7, —a coincide. Hence the double index 
7, a can be combined into a single one, a = (+, — ), with a = + for A f = B { and a = — 
for A I— B "J". Thereby the number of independent parameters is reduced. Eq. (|12"D then 
yields two coupled sets of self-consistent equations for G an and E an 

°? N(E) 
G an = / rfE: ■ ==jp ■ -; r (13a) 

— OO 



G an = - J dr e^ r ((^(r)C(r + ))r)a, (13b) 



where \x a = /i + ah s . Here (■ ■ -)t represents the statistical average with the local action as 

(0) T = -^- J Vif;VrOe A + +A - (14) 

Note that ( |13a|) merely expresses the fact that G an is the local element of the full Green 
function, i.e. is an explicit function of the self energy S an . Only eq. ( |13b| ) describes the 



actual dynamics of the system, determining the self-energy S an as a function of the input 
parameters T, /x, U and the disorder distribution. Eqs. (0) fully determine the one-particle 
properties of the Anderson-Hubbard model in d = 00. 

III. CORRELATION FUNCTIONS AND THERMODYNAMIC STABILITY 

To decide about the stability of a particular solution of ([i~3|) we have to evaluate averaged 
two-particle correlation functions or susceptibilities. Two such quantities are of particular 
interest in the Anderson-Hubbard model: the averaged staggered magnetic susceptibility 
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Xav 9 = Xaf an< i the averaged compressibility n av . The former decides about the instability 
of the paramagnetic phase w.r.t. AFLRO, while the latter not only contains the information 
about a possible instability w.r.t. phase separation (n av — » oo), but also provides a thermo- 
dynamic criterion for a solution to be insulating (n av — > 0). Quite generally susceptibilities 
can be obtained from the second derivative of the averaged grand potential Q av w.r.t. to 
some variable x as 



Y 1 92Q ™ 



XAFi % — h s 



E f X a G an (15) 



L dx 2 
The first derivative of Q av w.r.t. x yields 

1 on^ _ _^_ 

L ox 2 an 

where f£ = 1 and f^ s = a. Differentiating once more we obtain the desired susceptibilities 
(15). Note that these susceptibilities are evaluated at zero staggered field h s , i.e. h s is put 
to zero once the second derivative of fl av w.r.t. \i or h s is taken (in this case fi a = /i). Thus 
the susceptibilities are given by 

Xav =13 2^ fcXnn',n'nla'n' (16) 

an,a'n' 

where T is the local two-particle correlation function 



-pact 

n-\n' .nLnn 

1> 2 / 



/ dr x dr[dr 2 dT 2 exp {i(w ni Ti + uv^ - o; n ' 2 r 2 - uj m T 2 )} 



X [(^«(TiX(r 3 )^ W)^(^)>r>«, - ({Mn)ra(T2))T{M<)rM))T)ac] (17) 

The quantity 7^ n = 9(G r ~^ + E Q , n )/9x in (|T6|) measures the response of the averaged medium 
to an infinitesimal change of the field x. We will later see that 7^ n decides about the (in-) 
stability of a given phase. This dynamical response function is determined by an integral 
equation in frequency-space which does not explicitly depend on momentum. (Note that 
there are no convolutions in fc-space in the d = oo limit as is typical for a mean-field theory) . 
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This property does not imply, however, that the response function 7^ n is local, too. It only 
indicates that 7^ n is diagonal in the momentum q. A momentum dependence will enter by 
taking derivatives, e.g. of the local Green function, w.r.t. an external field with a particular 
g-dependence (q = in the case of the compressibility or the ferromagnetic susceptibility, 
and q — (it, . . . , ir) for the staggered susceptibility). Infrared divergencies show up in the 
spectrum of 7„(£). 

The integral equations for 7^ n are derived as follows. Using the effective propagator 

y m = G an "t" ^an = (1 + G anX-'an) / G am SUCh that 

lln = dGal/dx (18) 

we replace the self-energy in (13) by E an = Q~^ — G~^ n . Then (13) reads 

G an = (g a n(E)) E (19a) 

Gan = - f P dTe^ T ((MrWa(r + )) T )a V (19b) 



where g an {E) = (z an - E 2 / Z- an ) \ with z an = iu; n + fi + ah s - Q a £ + G a „, and 

(y(E)) E = / dEN(E)y(E). (20) 

J CO 

To obtain 7^ n we differentiate (|19|) w.r.t. x. This yields 



^an ~ I Ian ' r^l Jot I ^ an ' I ^ '—an ' r^2 J -a I Van [ZlB,) 



"an .«. \ *. / o. Ur 

G 'an = I 3 ' 1 E r nn',n'n7a'„' ( 21b ) 



where G' an = dG an /dx and £ an = (gl n {E)) E ,r] an = (E 2 g 2 an (E)) E / z 2 _ an . Eliminating G' an 
from (|2T|) we obtain the following matrix equation for 7+ n and 7^ n : 

/3- 1 £ (£<Wi^' + C^ Jt*v = <n (22) 

a'n' 

The quantities i?" a and v x an are a matrix and a vector, respectively, w.r.t. the sublattice-spin 
index a and may be written as 



13 



R„ = -(cietD n )- 1 D n T ri 

R-n " f X 



"J in 



(23a) 
(23b) 



where 



U-£~ 



D r 



cv- 



V-,n 



or: 



V Gl 



(24a) 



gTu/ 



Z,+,n T]+,n 



\ V-,n £-,n ) 



(24b) 



/l \ 



r 



(24c) 



1 o an 



= f- a vl an , eqs. (23)-(25) reduce to (Bf) in Ref. 



In the paramagnetic phase, where f*v 2 
[42]. 

To determine the boundary between the paramagnetic and the antiferromagnetic phase 
it is sufficient to calculate xaf only in the paramagnetic phase. By contrast, n av has to be 
evaluated in both phases since in both phases a metal-insulator transition (MIT) is possible, 
in principle. Indeed, in the paramagnetic phase of a random alloy with discrete disorder 
spectrum a MIT due to band-splitting may occur (at U — this is known to be the case, at 
least within CPA [[|0],[|T]), and even in the presence of AFLRO there are indications that both 
a metallic and an insulating phase exist |42j| . We will find that the most sensitive indicator 
for an incipient transition are the dynamic response functions 7^ n , x = [i,h s - especially their 
behavior close to the Fermi energy (i.e. their values for the lowest Matsubara frequencies), 
rather than \ AF and K av itself. The latter quantities are much less sensitive to changes in 
the ground state since they represent sums over the Matsubara frequencies. 



14 



IV. SPECIFICATION OF THE MODEL AND OF THE NUMERICAL 

PROCEDURE 

Eqs. (15)- ( pip form the basis for the numerical evaluation of n a v and Xaf- Before they 
can be solved quantitatively we have to specify the model parameters we use. The numerical 
calculations were performed with a semi-elliptic density of states (DOS) with total width 
2u>, i.e. 

N{E) = —(w 2 -E 2 ) 1 / 2 . (25) 

7rar 

This DOS is chosen because of its sharp algebraic band edges, resembling those typical for 
d — 3, and its simple analytic form; it is exact for a Bethe lattice in the limit Z —* oo. 

To study the influence of the disorder we investigate, and compare, three qualitatively 
different distributions of random potentials: 
a) Discrete, binary-random-alloy distribution 

, , 1./ An 1 / A 



t-binaryKti) = -<J(Cj - —) + -8\€i + —). (26a) 

The atomic potentials 6j = ±A/2 are chosen with equal probability to ensure an average 
band filling of n = 1 (only in this case, or n ~ 1, is AFLRO expected to occur at all). This 
distribution is important since it leads to a disorder-induced MIT due to band-splitting in 
the noninteracting system (an exact result in d = oo [38,39] ) which may compete with the 
interaction-induced Mott-Hubbard MIT. 
b) Continuous, semi-elliptic distribution 



*, semil iy 



VTA 2 



An2 



2 J - e ' 



1/2 



(26b) 



This is a much softer type of disorder. We found that a continuous, constant disorder 
distribution P cons t(ei) = A _1 ^(^- — |ej|), leads to essentially identical results, 
c) Percolation-type disorder 

PperM) = (1 - x)5(e l ) (26c) 
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By using P pe rc{u) one may simulate doping in a non-random system with impurities which 
do not hybridize with conduction electrons: with probability x an infinite energy barrier is 
created which prevents the electrons from visiting these sites. The three distributions allow 
us to test the universality of the magnetic behavior of the model obtained for different types 
of disorder. 

To be able to study the competition between magnetic order caused by the electronic 
interactions and the disordering effects caused by the random potential, respectively, we 
work with an average band filling n av = 1. Due to the symmetry of the distributions (|26a| ,b) 



we can fix the chemical potential at /i = U/2. This can even be done with ( |26c|) in which 
case the Green function has to be weighted with an additional factor of 1 — x. 

For the numerical evaluation of the functional integral (|14"D we employ the algorithm of 
Hirsch and Fye p8| . We discretize the time variable, i. e. j3 = ASt, with 0.25 < St < 1, 



and then extrapolate the quantities under investigation to St — > 0. To obtain a smooth 
imaginary-time Green function G(t) even for discrete values of r we use instead of G n (the 
Fourier transform of G(r) w.r.t. Matsubara frequencies u n ) the function 

G n = 5t/[1 - exp(Sr/G n )] (27a) 

G(t) = /T 1 J2 e iulnT G n (27b) 

n 

with (3 = ASt. In the continuum limit St — > the functions G n and G n coincide. The 
definition fl2"7|) preserves the relation G(t = 0) = n — 1 even in the case of discrete St 
and suppresses unphysical and undesirable oscillations in G(t) due to the discretization (see 
Fig. 1). In particular, by smoothing G(t) in this simple way large-scale Fourier transforma- 
tions (namely for every iteration) described in Ref. [^9j become unnecessary. 

Exact summations over spin variables in the discrete Hubbard-Stratonovich transforma- 
tion were used whenever possible, i. e. for A < 22. For A > 22 we used Monte-Carlo 
sampling. After typically 4-8 iterations an accuracy of 10~ 5 and 10~ 3 was reached in the 
exact summations and in the Monte-Carlo sampling with 10 4 sweeps per iteration, respec- 
tively. Close to the magnetic transition the convergence becomes significantly slower and 
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the Monte-Carlo sampling less and less efficient. This implies that in the immediate vicinity 
of the transition this method cannot be used to obtain accurate results. 

The integration over the continuous disorder distributions was performed using the 
Gauss-Legendre quadrature p5| . A discretization of the random energies in steps of 



Se/w ~ 0.2 gave an accuracy of ~ 10~ 5 for all observables. The required CPU time increases 
proportional to A 2 2 A and A 3 for the exact summation and the Monte-Carlo sampling, re- 
spectively ||52| . Although the random energies break the particle-hole symmetry we never 
encountered a minus-sign problem. 

Setting h — ks — 1 the only remaining physical dimension is that of an energy (U, E, A, T, 
etc.) or inverse energy (G, n av , XafiP-, St). Departing from our earlier convention we now 
choose the half-band width w as our energy unit since it does not depend on the limiting 
process Z — > oo. (This is in contrast to the scaling of the hopping amplitude t = at* /vZ, 
where a may be chosen at will). This convention agrees with that used by Kotliar and 



collaborators pl|j23|j30|j32]] . To be able to compare the results presented in this paper with 



our earlier ones [33|, all numerical values of quantities with dimension of energy (inverse 



energy) obtained in Ref. [|IJ must be divided (multiplied) by a factor of 2. To compare 



with the results of Jarrell's group P^J2^|J53[| and that of Georges and Krauth [^,[J]]], their 
numbers have to be divided (multiplied) by a factor of \/2. 

V. MAGNETIC PHASE DIAGRAM 

The competition between disorder and correlations in the Anderson-Hubbard model 
strongly affects its properties both in the ground state and at finite temperatures. Con- 
cerning the thermodynamics the influence of different kinds of disorder on the stability of 
AFLRO near half filling is of particular interest. To determine the instability of the paramag- 
netic phase w.r.t. the formation of AFLRO we evaluate eqs. (15) - (25) in the paramagnetic 
limit for x = h s , in which case the dependence on the index a drops out; we set S an = E n , 
etc. and 7^ = j n . Note that the matrix equation (^) separates into two identical scalar 
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equations in this case. As already mentioned below (|17|) it is not the averaged susceptibility 
itself which is of primary interest in the investigation of the stability of the paramagnetic 
phase but the dynamic response function j n . The averaged susceptibility diverges (becomes 
negative) if and only if the real part of the response function 7„ diverges for at least one fre- 
quency. This is characteristic for the dynamics of a quantum system of interacting particles 
possessing infinitely many coupled internal degrees of freedom labelled by the Matsubara 
frequencies. Each frequency corresponds to one mode in the quantum mechanical system 
which is described by complex variables. Note, however, that the Matsubara frequencies do 
not index the actual independent modes, since T nn / >n / n is not diagonal in frequency. Only 
the eigenvalues of the integral equation (^2|) represent the independent (eigen-) modes of the 
interacting quantum system. We can ascribe a critical (Neel) temperature T/v,|„| to each 
eigenmode n whereby T Nt \ n \ > Tjv,| n '| if \ n \ < l^'l- in Fig- 2 we show the function 7 n (r), 
being a linear combination of eigenmodes for the lowest Matsubara frequencies which lie 
very close to the diverging independent soft mode. The functions 7„(t) diverge due to the 
contribution of the eigenmode with the highest T N \ n \. The highest critical temperature T N \ n \ 
is the thermodynamic critical temperature T^. 

Once the function 7„ and the local vertex function T nn >y n are known one may calculate 
the averaged susceptibility. Above T/v this susceptiblity must obey the Curie- Weiss law, 
provided there is such a T/v > 0. In Fig. 3 three characteristically different temperature 
dependencies of Xaf are shown for one value of (binary alloy) disorder, A = 2, at different 
interactions strengths. If U/A is sufficiently large the Curie- Weiss law with a finite T/v is 
obeyed as in the case without disorder. For values of U ~ A the low-temperature behavior 
begins to be determined by the scattering off the frozen random configurations and deviations 
from the Curie- Weiss-law become apparent. For (//A < 1 we observe a minimum in Xaf 
(i.e. a maximum in xaf itself) which separates the temperature-dominated regime from the 
disorder-dominated regime. In the latter case the long-range correlations are continually 
suppressed due to impurity scattering. The phase diagram in the T — ?7-plane calculated 
from the zeros of Xaf i s plotted in Fig. 4a for different values of the binary alloy disorder. We 



can distinguish two different regimes, a) U < U c g 2.5 : the disorder gradually suppressed 
the long-range order and thus reduces the critical temperature. This is exactly what one 
expects from scattering off frozen random configurations; b) U > U c : here the situation is 
strikingly different since T/v-curves of constant disorder start to cross. This means that a 
small amount of disorder supports the formation of AFLRO, i.e. the critical temperature 
increases with disorder as shown in Fig. 5. A maximal critical temperature is reached at 
some value A C (U) beyond which a further increase of disorder causes the critical temperature 
to decrease monotonically to zero. This effect is particularly pronounced at strong coupling. 
It is also observed, but less pronounced, in the case of the continuous disorder distribution 



(Fig. 4b). In the limit U ^> A, t the enhancement of T/v may be explained as follows [54 
The virtual hopping of an electron with spin a from a given site A with local energy e A to 
a neighboring site B with energy e B occupied by a (— a) electron leads to an energy gain 
J\ = —t 2 /[U — {ca — e B )\ and to J 2 = —t 2 /[U + (e A — €b)] for the reverse process. The 
effective spin coupling J(e A , e B ) is given by the sum of these energies, J(e A , e B ) = J\ + J2, 
which for a bounded disorder distribution, — ■§■ < £a,b < f- , with U 3> A (strong coupling) 
implies 

2t 2 



J(^A,e B ) 



1 + 



fe A -e B ^ 2 



(28) 



V u 

Assuming that even in the presence of disorder the Neel-temperature T^(U,A) 
oc< J(e A , e B ) >av, we find 

r^0) } = J deA J deB J(^,e B )P(e A )P(e B ) = 1 + -^[p 2 - p\] (29a) 

for arbitrary disorder distribution P(e). p\ and p<i are the first and second moment of P(e), 
where p\ = J dee P(e). For a symmetric, bounded distribution (pi = 0) one therefore finds 

T ^ A) -1 + A^ 2 (29b) 



T N (U,0) \U. 

with A = 2p 2 /A 2 , i.e. the disorder is indeed found to increase T/v irrespective of the type 
of disorder. This effect is the more pronounced the more structure P(e) has at |e| ~ A 



2 ■ 
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i.e. is larger for the binary alloy (A = ~) than for the semi-elliptic (A = |) distribution. 
This quadratic increase of T/v with A is indeed found numerically at large U (Fig. 5). 
Apparently the mechanism that enhances T/v is effective already at ^/-values as small as 
U ~ 2.5. It is interesting to note that for given strength of disorder A the [/-value where 
AFLRO begins to set in, U m i n , given by Tjv(C/ m j n , A) = 0, tracks almost perfectly with A, 
i.e. U min (A) ~ A. More generally, for binary alloy disorder the T/v-curves are found to obey 
the phenomenological scaling law 



T N (U,A)~T N [U 



,0 



2 \U 



(30) 



for arbitrary U (Fig. 6). For U ^> A (31) reduces to (29b). 

It is well-known that AFLRO may be destroyed by doping the (non-random) system with 
holes, such that n av < 1. This is a very interesting effect both from an experimental and 
theoretical point of view. The effect of adding holes may be approximated by introducing 
percolation-type disorder, ( |26c| ), where the random site-energy is zero with probability 1 — x 
and infinite with probability x. Since in d = oo the critical value for percolation is x — 1 



55| AFLRO will persist up to x — 1. The boundaries between the paramagnetic and the 
antiferromagnetic phases for this type of disorder are shown in Fig. 7. At weak coupling 
the critical temperature remains constant over a wide range of x and then rapidly falls to 
zero. At strong coupling, however, the critical temperature monotonically decreases and 
approaches the linear behavior T/v(x) = T/v(0)(l — x). The latter dependence of T/v which 
is very different from the one observed in the previous two cases is an exact result in the 
case of the disordered Falicov-Kimball model |39| . 



The magnetic correlations are mediated by the local magnetic moments. Their static 
average, m av , is defined by m 2 av = I/ -1 Ei(((^iT _ ^|) 2 )t )av = 1 - 2d av . Here d av = 
L~ l dQ av /dU is the average double occupancy of lattice sites. For the binary and the semi- 
elliptic disorder distributions discussed above m 2 av is shown as a function of U in Fig. 8a. 
The disorder is seen to have two main effects which are independent of the specific disorder 
distribution: 1) at a fixed value of U an increase of the disorder reduces the moments. 
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2) For fixed disorder strength A an increase of U leads to an increase of the moments, with 
saturation starting at U ~ A. Both features are easily explained in terms of the effect 
the disorder and the on-site repulsion, respectively, have on the average double occupancy 
dav For (/ -C A the local repulsion is weak while the spatial fluctuations of the atomic 
potentials are strong, such that d av is at its maximum value (d av ~ 1/2) and m av is small. 
As U increases the particles are forced to separate and d av decreases, i.e. m av increases, 
too. The local magnetic moments do not show a critical behavior close to the transition 
temperature T/v. In the case of the percolation-type disorder (Fig. 8b) the density of local 
moments on the reduced lattice with (l — x)L sites increases with concentration x. In other 
words, in the presence of disorder, the saturation of m 2 v sets in at smaller values of U than 
in the pure system. The influence of the percolation-type disorder can be explained by the 
reduction of the kinetic energy by a factor of 1 — x which is due to the reduced average 
number of nearest neighbors. For large U, in the Heisenberg limit, the factor t 2 in the 
antiferromagnetic coupling, J ~ t 2 /U, has to be replaced by t 2 (l — x) leading to a linear 
decrease of T/v- For small values of U, the influence of the repulsion increases due to the 
reduction of the kinetic energy. Therefore, the double occupancy is suppressed, i.e. m 2 av is 
enhanced by the percolation disorder. The enhancement of m 2 v leads to the stabilization of 
AFLRO. This effect is definitely not described within the Hartree-Fock approximation; the 
latter only leads to an exponential suppression of T/v if U or (1 — x) becomes small. 

In Fig. 9 the temperature dependence of the order parameter (the averaged staggered 
magnetization Maf) is shown. In the vicinity of the transition point, T ~ T/v, it is well 
represented by a mean-field- type dependence Maf(T) oc (T n — T) 1 / 2 indicated by the dashed 
line. We observe that the extrapolation of the Monte-Carlo data to Maf{Tn) = using 
this law leads to the same critical temperature as that obtained from the divergence of Xafi 
i.e. Xaf(Tn) — 0, assuming a Curie- Weiss law for T ~ T N . Deviations from the square-root 
behavior of Maf{T) set in at lower temperatures. 
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VI. METAL-INSULATOR TRANSITIONS 

In the preceding section we found that at weak coupling the disorder suppresses the 
AFLRO. In the ground state the system is then a disordered paramagnet. Whether it is 
metallic or not depends on the type and strength of the disorder. In d = oo, where Anderson 
localization does not occur [38,39] a paramagnetic insulator only forms if the spectrum of the 
disorder distribution is multiple connected, such as in the case of the binary-alloy disorder. 
For this type of randomness a MIT due to band splitting is expected to occur at some value 
A > 1. In the non-interacting case, U = 0, the exact result in d = oo (obtained by CPA 

m)isA=i. 

For U = there are two equivalent criteria to decide on whether the ground state in 
d = oo is metallic or insulating. They correspond to (i) a spectral definition of an insulator, 
based on the disappearance of the DOS at the Fermi level (assuming that only extended 
electrons are present), and (ii) a thermodynamic definition, employing the disappearance 
of the compressibility [|5(J]. However, in interacting systems these two criteria need not 
coincide. Hence we will investigate both. In Fig. 10a the averaged compressibility n av of 
interacting electrons in the presence of binary-alloy disorder is shown as a function of U 
for (3 = 16. Hence K av is calculated across the line T = 0.0625 in Fig. 4a, i.e. within 
the paramagnetic and the antiferromagnetic phase, respectively (the transition points are 
indicated by arrows). For a given A ~ 1 the curves display a common behavior as a function 
of U: K av is (exponentially) small at U = 0, then increases and, at U ~ A, approaches a 
maximum, beyond which it becomes (exponentially) small again for [/> A. This behavior 
has a clear physical interpretation: For U <C A, with A ~ 1, the double occupancy of 
lattice sites is at its maximum value (d av ~ 5), at the same time the kinetic energy cannot 
delocalize these states. Hence the system is insulating. At T = and U = this is an 
exact property for A > 1 (split-band limit) in d = 00. As U increases the on-site Coulomb 
repulsion forces the particles to separate from each other. As a consequence the particles 
become less localized. We then expect that, at some critical value U^' 11 ' 1 , a macroscopic 
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fraction of the electrons becomes extended and the system starts to be metallic. This real- 
space picture has its analogy in A;-space (or w.r.t. energy): as the interaction increases the 
previously separated bands change their shapes, i.e. the interaction-induced energy exchange 
between particles leads to a transfer of states into the energy gap where the Fermi energy 
is located. At the critical value U^ 11,1 the (still algebraic) band edges reach the Fermi level, 
producing a finite DOS there as well as a finite overall compressibility. (We cannot rule out 
that U^ 11 ' 1 = 0, i.e. that the DOS at the Fermi level is finite even at arbitrarily small U; 
however, this would require the bands to acquire exponential tails - a feature which cannot 
be observed numerically). Hence, in contrast to an interacting system without disorder the 
Coulomb interaction in a disordered system is able to improve the metallicity of the system 
and may even turn an insulator into a metal. Apparently the interaction-induced energy 
exchange smooths the energy spectrum of the disordered system, thus leading to an easier 
transfer of energy. 

For U > U^ 11 ' 1 the system is then expected to be a metal without Fermi liquid properties 
(since JmE(w = 0) 7^ 0), at least in d = 00. As U is further increased the DOS and K av 
increase, reaching a maximum at U ~ A. For even larger U the effect of pushing states into 
the gap is reversed and K av decreases again. Hence at some critical value U^ 41,2 a second 
transition occurs, back into an insulating state, where doubly occupied sites are almost 
completely suppressed, whereby the mobility is obstructed by the repulsive interaction. 

In Fig. 10b K av for the continuous, semi-elliptic disorder distribution is shown. Here 
K av behaves qualitatively as in the non-random case, i.e. n av decreases monotonically since 
the DOS does not split at U — 0, T — 0. This reduction of K av occurs in two steps: for 
U ~ A K av decreases almost linearly with U, corresponding to the linear suppression of 
doubly occupied sites d av , or m 2 v (see Fig. 8a); in this region the system is metallic but not 
a Fermi liquid. Then, for U ~ A, when the number of doubly occupied sites is almost zero 
so that the local moments are almost saturated, K av approaches zero (exponentially) slowly. 
Strictly speaking, at finite temperatures one has n av > for U < 00. 

Next we discuss the relation of U^ 41,1 and U^ 41,2 to the critical interaction strength 
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U^ F where the paramagnet becomes unstable w.r.t. AFLRO. In Fig. 10a,b the position 
of U^ F is indicated by an arrow. Since the local moments vanish in the (binary-alloy)- 
disorder-induced insulating phase, the insulator-to-metal transition at U^ 1,1 occurs before 
the magnetic order sets in, i.e. U^f 1,1 < U^ F . This raises the question of what happens as one 
goes through the magnetic transition into the ordered phase: will the compressibility jump 
to zero discontinuously or will there be an antiferromagnetic metal, i.e. does the interaction- 
induced MIT coincide with the magnetic transition or not? To answer this question the 
behavior of K av near the transition point has to be investigated. In Fig. 11a the change of K av 
due to the onset of AFLRO in the case of binary-alloy disorder is shown. The results for the 
paramagnetic phase are compared with those for the ordered phase. At A = the AFLRO 
is seen to suppress K av drastically. By contrast, at finite disorder (A = 2), this difference 
almost vanishes. Apparently the disorder stabilizes the metallic state close to U^. A similar 
behavior is observed for the continuous disorder (Fig. lib). These results suggest that in 
the vicinity of Uf F the system is an antiferromagnetic metal. However, MC-techniques are 
not able to decide whether this is true even at T = 0. We attempted to extrapolate our 
finite-temperature data to T = but could not find a simple, accurate extrapolation law. In 
Fig. 12 the averaged compressibility for the system with binary-alloy disorder, calculated at 
T = 1/16 and 1/40, is shown. Although at the lower temperature K av is lower and the slope 
is slightly steeper at the transition point, the critical point (indicated by an arrow) seems 
to be well inside the metallic phase. In particular, the tail behavior for U > U^ does not 
show any significant change. 

To investigate whether, and how, an insulator evolves it is instructive to plot the dy- 
namical response function 7 an , or the quantities 

K « := o H dG an /d/j, (31a) 

with 

Kav =l3~ 1 ^2^n, (31b) 

n 
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as a function of Matsubara frequency u n = irT(2n+l). In Fig. 13a the real part of K n is shown 
for a system without disorder, A = 0, for U = 1.75 and T = 1/64 (the system is then well 
inside the antiferromagnetic region, see Fig. 4a). The behavior of ReK n close to u n = is seen 
to be very different for the antiferromagnetic and the (hypothetical) paramagnetic solution. 
For the paramagnetic solution, which is metallic, Ren n is a monotonically decreasing function 
of \oo n \. By contrast, in the antiferromagnetic (i.e. insulating) phase this is only so for 
\oo n \ ~ 1, while for uo n close to zero i?e«; n becomes strongly negative. In the former case the 
sum over n n , i.e. K av , is then clearly positive, while in the latter case K av becomes very small 
due to the negative contributions close to u n = 0. This behavior at A = is now contrasted 
with that in the presence of binary-alloy disorder (A = 2, see Fig. 13b) at U — 2.65 for 
two different temperatures: T = 1/20 (close to the magnetic transition) and T = 1/64 
(well inside the ordered phase; see Fig. 4a). At T = 1/64 the range of u n values for which 
Ren n < is now even narrower the in the case without disorder. This shows that a definite 
answer to the question of whether the antiferromagnetic phase close to U^ F is insulating 
(n av = 0) or not at T = can only be obtained from the behavior of Ren n at very small 
u n , i.e. from K av at very low temperatures (T < 10~ 2 ). Within the MC-approach used here 
such low temperatures are not attainable. 

The situation is similar in the case of —ImG n , the imaginary part of the one-particle 
Green function G n = | J2 a G an , which yields a spectral condition for an insulator (for n = 
and in the limit T — > the function —IitlGqiv coincides with the DOS at the Fermi level). If 
ImGo — ► for T — >• the ground state of the system is insulating. In Fig. 14 the influence of 
the onset of AFLRO on —IitlGq is shown. Although in the disordered system the difference 
between the magnetic and non-magnetic solution is now slightly stronger than in the case 
of K av (see Fig. 13a) the magnetic transition occurs where ImGo is large, i.e. well inside the 
metallic phase. In Fig. 15 the behavior of —ImG n at the few lowest values of u n > is 
plotted for the same values of U and (3 as for Ren n in Fig. 13. In the case without disorder 
(Fig. 15a) this behavior is qualitatively similar to that of Ren n (Fig. 13a): the paramagnetic 
solution is clearly metallic while the antiferromagnetic phase is insulating since IitlGq — ► 0. 

25 



In the presence of disorder (Fig. 15b), however, the downturn of —ImG n in the phase with 
AFLRO (/3 = 64) sets in only very close to \u n \ ~ i.e. at very low temperatures. If close 
to the transition the antiferromagnetic phase were really insulating the decrease of —ImG n 
will have to be very rapid. In Fig. 16 we plotted our Monte Carlo data for —ImGo in the 
case of binary-alloy disorder with A = 2 as a function of temperature for several values of 
U. Both at very small U (i.e. U < 1, where the disorder- induced gap of the DOS dominates 
the behavior) and at large U (Heisenberg-limit) the temperature dependence may be clearly 
extrapolated to —ImGo = for T = 0, i.e. the ground state is insulating. However, for 
U ~ U^ F (A,T = 0), with U^ F (2,0) ~ 2 — 2.5, these results cannot be safely extrapolated. 
Hence the question whether in the disordered system close to the magnetic transition the 
antiferromagnetic phase is metallic or not remains open. As in the case of n av one would 
have to go to temperature T ~ 10~ 2 to be able to decide whether the ground state of the 
disordered system with AFLRO is metallic or insulating. 

VII. CONCLUSIONS 

In this paper we presented a detailed, quantitative study of the physical effects caused by 
the simultaneous presence of interactions and randomness in a system of lattice electrons. 
To this end we investigated the Anderson-Hubbard model with diagonal disorder at half 
filling in the limit of infinite spatial dimensions, i.e. within a dynamical mean-field theory, 
for three different disorder distributions. Numerical results were obtained by employing 
Quantum Monte Carlo techniques that provide an explicit finite-temperature solution of the 
model in d = oo. No further approximation were used. 

To construct the thermodynamic phase diagram we derived and evaluated the appro- 
priate averaged two-particle correlation function, i.e. a dynamical response function, whose 
poles determine the magnetic instabilities of the disordered, interacting system. Only this 
function - and not the averaged staggered susceptibility Xaf itself, which is only a weighted 
sum of the response function over the (Matsubara) frequencies - fulfills a closed equation 
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that determines the two-particle spectrum. The value of the response function at the lowest 
frequency is the most sensitive indicator for an instability of the system. In the temperature 
dependence of Xaf two distinct disorder regimes are observed: (i) for weak disorder the 
Curie- Weiss law holds, while (ii) at strong disorder Xaf acquires a maximum at a tempera- 
ture below which a crossover from the temperature- to the disorder- dominated regime takes 
place. 

We demonstrated that at low temperatures and sufficiently strong interaction there al- 
ways exists a phase with antiferromagnetic long-range order (AFLRO). Furthermore we 
discovered a new strong-coupling anomaly, namely that the Neel-temperature T/v is not al- 
ways a monotonously decreasing function of disorder. Indeed, at strong coupling and not 
too large disorder T/v is always found to be an increasing function of disorder, i.e. disorder 
favors the formation of AFLRO in this regime. This implies the existence of an unexpected 
disorder-induced transition to a phase with AFLRO. Under the assumption that T/v is pro- 
portional to the effective exchange coupling between spins even in the disordered system we 
proved that for diagonal disorder the anomalous behavior is generic, i.e. is independent of 
the type of disorder distribution. It is a consequence of the fact that for diagonal disorder 
the difference between the local energies of neighboring sites becomes smaller on average, 
thus leading to stronger effective exchange coupling. 

We then studied the existence of metal-insulator transitions in the Anderson-Hubbard 
model. Although in d = oo Anderson localization does not take place the presence of 
disorder may well have other strong effects. In particular, binary-alloy disorder is able 
to cause band-splitting (thereby resembling the effect of genuine interactions) and hence 
may induce a metal-insulator transition all by itself. Special attention was given to the 
question whether or not the disorder allows for the stabilization of an antiferromagnetic (AF) 
metal. To this end the average compressibility n av was evaluated both in the paramagnetic 
(P) and AF phase. Contrary to our expectation the presence of disorder was found to 
enhance the metallicity of the AF-phase close to the P-AF transition. This enhancement 
strongly suggests the existence of an AF metal in the low-temperature phase of the Anderson- 
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Hubbard mode. To investigate whether this AF phase persists to be metallic down to T = 0, 
at least in d = oo, we studied the frequency components of the one-particle Green function 
and of K av , respectively, down to T = 1/64, our lowest temperature. However, a reliable 
answer to this question can only be found at still lower temperatures, which at present 
are beyond the reach of the finite-temperature Monte Carlo techniques used here. In any 
case, the transition scenario involving metal-insulator and P-AF transitions obtained for the 
Anderson-Hubbard model is remarkably rich. For alloy-type disorder with A > 1 (split-band 
limit for U = 0) at T = an increase of the interaction U from zero will probably first lead 
to a transition from a paramagnetic insulator to a paramagnetic metal at U^ 11 ' 1 , then, at 
U^ F , to an antiferromagnetic metal and finally, at U^ 1,2 , to an antiferromagnetic insulator. 
No compelling evidence was found that U^ F and Ujf 1 ' 2 coincide. In the case of a continuous 
disorder distribution one has Uff 1,1 = since band splitting never occurs. 

The above findings prove that the interplay between electronic interactions and scattering 
from disorder leads to interesting, and even novel, physical effects. In particular, the strong 
coupling anomaly discovered here call for an experimental verification. 
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Figure Captions 

Fig.l Local Green function G(r) for U = and /i = (a), fi = 0.5 (b) obtained by usual 
Fourier transformation (dotted line) and by redefinition ace. to eq.(13) (solid line) 

Fig. 2 Dynamical antiferromagnetic response function 7„ vs. temperature T. Binary alloy 
with disorder strength A = 1 and U = 2. Here and in the following figures lines are 
usually guides for the eye, and error-bars are roughly of the size of the symbols unless 
shown explicitly. 

Fig. 3 Inverse averaged antiferromagnetic susceptibility Xaf vs - ^ f° r the binary alloy with 
A = 2 and several values of U. 

Fig. 4a T — [/-phase diagram for the binary alloy with A = 0,1,2,4 obtained from the 
zeroes of Xaf ( see Fig. 4). The AF-phase is stable below the curves. The dotted 
lines at T = depict the regimes where the Curie law would give negative transition 
temperatures. Below the crosses Xaf nas no zeroes but a minimum and an AF-phase 
can no longer be expected. 

Fig. 4b T — [/-phase diagram for the semi-elliptic distribution of the random energies with 
width A = 0,2, 4, 6. 

Fig. 5 T vs. A for the binary-alloy distribution at U — 3, 4 and 5.5. Dashed lines: quadratic 
increase of T/v according to eq. (30b) with A = 1/2; dotted lines are guides to the eye 
only. 

Fig. 6 Scaling plot of T/v for the binary-alloy distribution according to eq. (30); symbols as 
in Fig. 4a. 

Fig. 7 T — x-phase diagram for the percolation- type disorder for several values of U. The 
dashed line depicts the Hartree-Fock (HF) result for the smallest value of U — 0.65. 
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Fig. 8a Averaged quadratic local moment m av vs. U at inverse temperature /3 = 16. Without 
disorder (full circle); binary alloy with A = 2 (full square); semi-elliptic distribution 
with A = 2 (open circle) and A = 4 (open square). 

Fig. 8b Averaged quadratic local moment m 2 av /(l — x) (normalized to the concentration of 
sites with random energy zero) vs. x for the percolation-type disorder at /3 — 16 for 
several values of U . 

Fig. 9 Staggered magnetization Maf and inverse averaged susceptibility Xaf vs - T f° r the 
binary alloy with A = 1, U = 2, St = 2/3. Dashed line: square-root fit of the last two 
points of Maf below the transition; dotted line: linear fit of Xaf- The arrow indicates 
the extrapolated Neel temperature. 

Fig. 10 Averaged compressibility K av vs. U at f3 — 16. a) binary alloy, A = 0,1,2,4; b) 
semi-elliptic distribution with A = 0, 2 and 4. Arrows indicate the transition to the 
antiferromagnetic state. 

Fig. 11 Averaged compressibility K av vs. U at (3 — 16 in the paramagnetic phase (dashed 
lines) and in the antiferromagnetic phase (dotted lines), a) binary alloy, b) semi-elliptic 
distribution. 

Fig. 12 Averaged compressibility K av vs. U for the binary alloy at /3 — 16 and 40. Arrows 
indicate the transition to the antiferromagnetic state. 

Fig. 13 Real part of K n vs. uo n , a) for A = 0, U — 1.75 and (3 = 64 in the paramagnetic 
phase (P) and the antiferromagnetic phase (AF); b) for the binary alloy with A = 2, 
U = 2.65 in the AF-phase at (3 = 40 and 64. 

Fig. 14 Imaginary part of the one-particle Green function G n at the lowest Matsubara fre- 
quency c^o = kT vs. U for a system without disorder (A = 0) and with binary-alloy 
disorder (A = 2) at (3 — 16; P-Phase (dashed lines), AF-phase (dotted lines). 
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Fig. 15 Imaginary part of G n vs. uj n in the P- and AF-phase; a) no disorder, U=1.75, f3 = 64; 
b) binary alloy with A = 2, U = 2.65, (3 = 20, 64 (the two curves at (3 = 20 in the 
middle are undistinguishable). 

Fig. 16 Imaginary part of Go vs. T for the binary alloy with A = 2 for several values of U. 
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